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Abstract 

We investigate the possibility of producing neutrinoless double beta decay without having 
an electron neutrino with a mass in the vicinity of 1 eV. We do so by having a much 
lighter electron neutrino mix with a much heavier (m > 1 GeV) sterile neutrino. We 
study the constraints on the masses and mixings of such heavy sterile neutrinos from 
existing laboratory, astrophysical and cosmological information, and discuss the properties 
it would require in order to produce a detectable signal in current searches for neutrinoless 
double beta decay. 



1. Introduction 



Our ultimate understanding of the universe relies on two crucial ingredients: the na- 
ture of the constituents (or building blocks), and the nature of the forces through which 
they interact. Experience to date indicates that spin-half fermions are the ultimate con- 
stituents, whereas forces arise from the exchange of gauge bosons associated with local 
symmetries, and Higgs bosons needed to break these symmetries. The highly successful 
Standard Model (SM) is based on 45 chiral fermions (15 for each generation consisting of 
twelve coloured quark states and three leptonic ones). Of these, the neutrinos have the 
unique property of being electrically neutral. 

A necessary part of the exploration of physics beyond the SM is the study of new 
chiral fermions and their properties. Apart from sheer curiosity regarding their existence, 
there are often very good physical motivations for postulating such particles. For instance, 
it is by now well known that if neutrinos have a mass, the smallness of this mass is easily 
understood if there is an additional heavy, isosinglet Majorana neutrino which mixes with 
the known ones via the so-called see-saw mechanism [1]. Similarly, there also exist quite 
good reasons which would motivate the existence of new charged fermions not present in 
the standard model. An important area of investigation in particle physics now is the study 
of the constraints that may be inferred on the properties of such new fermions using the 
existing data, as well as the identification of new experiments which can further constrain 
these properties. 

Since charged fermions cannot be singlets under the SM gauge group, the LEP data 
on the Z-width already imposes very stringent limits on their masses, i.e. m F > 45 GeV. 
However, no such direct constraints need apply for SM-singlet particles, unless they mix 
strongly with the known neutrinos. There are, however, many indirect limits on such 
particles [2], [3], [4], [5], [6], [7] and our goal in this paper is to explore these in the 
context of a specific model and indicate how future neutrino less double beta decay (/3/3q^) 
experiments can probe the existence of these particles in interesting ranges of masses and 
mixing angles. 

One motivation for examining the particular type of models we consider here is that 
they contradict the maxim 1 that the observation of neutrinoless double-beta decay would 
demonstrate the existence of a mass for the electron neutrino in the range of 1 eV [10]. In 

1 Alternative scenarios have also been explored for having observable f3(3 0l , decay without requiring 
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the models we study, all of the neutrinos are either much lighter than, or much heavier than 
the eV mass region. Neutrinoless double beta decay proceeds in these models through the 
virtual exchange of the heavier (i.e. GeV-scale or higher) neutrinos. We are able to find 
phenomenologically acceptable masses and mixings for such a model that are consistent 
with an observable signal for double beta decay. 

2. The Model 

We restrict ourselves for simplicity to supplementing the SM with two sterile neutrinos, 
N±, which mix only with the SM neutrinos of the first generation. The reasons for adding 
two sterile neutrinos rather than one are twofold: (a) the case of one extra sterile neutrino 
falls into a subclass of models which use the see-saw mechanism, and which have been 
discussed elsewhere [3], [4]; and (b) in models with one sterile neutrino, in the limit m Ue — > 
0, the sterile neutrino completely decouples and becomes invisible, whereas in models with 
two (or more) sterile neutrinos, the situation completely changes, and the sterile neutrinos 
can mix appreciably even in the limit m v< _ — > 0, and so can be potentially visible in many 
processes [5]. 

The interaction lagrangian involving N± can be written as follows: 

C = Co + d, (1) 

where C conserves lepton number and C\ violates it. (We assign L(N±) = ±1 to the 
left-handed parts of N±.) The most general renormalizable couplings and masses are: 

C = -mN +lL N_ - A e (L e7i iV_) H + h.c, (2) 

and 

d = N +lL N + - ^ N- lL N_ - A e (L elL N+) H + h.c, (3) 

where H is the usual SM Higgs doublet, and L e = (" e e ) is the first-generation lepton 
doublet. After the SM gauge symmetry breaking — with (H) = v = 174 GeV — one gets 

m„ e to be in the eV range [8]. Our approach here is more similar to, and updates, the framework of 
Ref. [9]. 
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the following most general mass term relating the left-handed states u e , N + and iV_ 



r = 
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~2 
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N + 









X e v X e v 



X e v 
X e v 



m 



m 



(4) 



We assume for simplicity, in writing eq. (4), that the elements of the mass matrix are 
real. This mass matrix has three nonzero eigenvalues, one of which is the physical electron 
neutrino, v' e . Since the direct laboratory bounds on the mass of a dominantly electron 
neutrino are quite low, we must ask under what circumstances a massless neutrino can 
emerge from this matrix. There are two possible regimes: 

(1) A e = fj, + = 0. The mass eigenstates can in this case be written as: 







(1! 
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SiC 2 CiC 2 — s 2 
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(5) 



where Cj (sj) denotes cos#i (sin^) while 



tan 9\ = 



m 



and tan 26 2 = 2 



\/ m 2 + X e 2 v 2 



(6) 



Finally, the factors of L V in the last row of eq. (5) come from the chiral rotation that 
is required to ensure that all of the entries in the final mass matrix are positive. The 
corresponding masses are: 



M v , = 0, 



and M N > = - 



(7) 



Notice that in the limit /j- — > 0, lepton number is exactly conserved, and iV4 forms 
a degenerate Dirac pair. In this case all lepton-number violating processes, such as 
I3[3q v decay, are forbidden. 

(2) X e = /j- = 0. This alternative is identical to case (1) above, but with the two sterile 
states N± interchanged. The physical implications of this case are therefore identical 
to case (1), and we need not further pursue this alternative separately. We henceforth 
exclusively focus on mass- matrix parameters that are in the vicinity of case (1). 
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For later purposes it is useful to define a dimensionless parameter, e, which measures 
the strength of the lepton-number violation in the mass matrix: 



e = = (8) 
V m 2 + A e v 2 

The ratio of the nonzero mass eigenvalues have a simple expression in terms of this pa- 
rameter: 

M N > Jl + 4/ e 2 + 1 

— — — — = v 7 = cot 2 9 2 . (9) 

M nV + 4/e 2 - 1 V ' 

In the lepton-conserving limit, e< 1, M N >_ ~ M^/ as would be expected for pseudo-Dirac 
particles. In this case the mass eigenstates, N'±, have opposite CP properties. For e> 1 
we instead have M N _/M N+ ~ e 2 , so that in this limit N'_ is much heavier than N' + . 

We could now turn on small non- vanishing values for A e and n+, so that the electron 
neutrino acquires a small nonzero mass. However, we are interested in the situation for 
which M Ve <C 1 eV and yet for which there are nevertheless potentially observable con- 
tributions to (3(3q u decay. We therefore keep A e and n+ negligibly small in what follows, 
although we return to the naturalness of this choice in the following section. We envision 
the possibility that the N'± masses span a wide range of possible values, starting from the 
MeV up to the GeV range. Clearly, the properties of such particles are highly constrained 
by known laboratory, astrophysical and cosmological information. We discuss the ranges 
of masses and mixings that are not already ruled out by these constraints and see if /3/3 0u 
decay can be observable in the allowed range. 

3. Contributions to f3/3 0l/ Decay 

As may be seen from eq. (5), the heavy sterile neutrinos acquire charged-current 
weak interactions through their mixing with the electroweak eigenstate v e , resulting in 
Kobayashi-Maskawa-type mixing angles 

Uev> = Ci, U eN > = SiC 2 , U eN > = -IS\S2- (10) 

e + — 

Provided that \i- ^ 0, so that lepton number is broken, neutrinoless double beta decay 
arises in this model via the exchange of the three neutrino states. The differential decay 
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rate for this decay between two + nuclei can be written in the following simple form [11]: 




cos 



|W| 2 8(Q -e 1 - e 2 ) [e lPl F(e 1 )}[e 2 p 2 F(e 2 )]. 



(11) 



ds\de 2 



2tt 3 



where G F is Fermi's constant; 9 C is the Cabbibo angle which governs the strength of the 
hadronic charged-current; Si and pi are the energy and momentum of each of the final 
two electrons; Q = M(Z,A) — M(Z + 2, A) — 2m e is their endpoint energy — typically 
several MeV; and F(e) is the Fermi function which describes the distortion of the electron 
spectrum due to the nuclear charge. Z and A represent the charge and mass number of the 
initial nucleus. For our present purposes it is convenient to work with analytic expressions 
for the total decay rate, which we can obtain if we make some simplifying assumptions, 
which are sufficiently accurate for the estimates in this paper. In performing the phase- 
space integrals we therefore neglect: (i) the electron mass, and (ii) the Coulomb-distortion 
factor, F{e). 

The dependence on the neutrino masses, in eq. (11), lies in the quantity W, which is 
given explicitly by: 



The sum here is over all three neutrino species, w = w(po, |p|) represents a particular 
Lorentz-invariant combination of form factors describing the nuclear matrix element of 
the two hadronic charged currents. All of the theoretical uncertainty in the decay rate 
enters with the estimating of w within a model of the nucleus. The connection between W 
as defined here and the usual estimates [12], [13], [14], based on an independent-nucleon 
model of the nucleus within the closure approximation, is given by [11]: 



w = Y^f± (N'(Z + 2,A)\J2t+t+ h{v nm - mi ) (gl - gla n ■ a m ) \N{Z,A)), (13) 



where m and n run over the labels of the nucleons in the nucleus; r+ and a m are isospin 
and spin matrices for the mth nucleon; and g v and g A are the vector and axial charged- 
current couplings of the nucleon. The function h(r mn ;mi) of the internucleon separation, 
r mn , is the neutrino potential, which is defined by the following integral: 




(12) 



mn 





(14) 
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where [i = \ [M(Z, A) + M(Z + 2, A)] is the mean excitation energy of the nucleus. 

For the present purposes, however, we need not use the detailed matrix-element ma- 
chinery, as quite good analytic results can be obtained by making the following simplifying 
approximation for the functional form for w. We parameterize w by representing it as a 
step functions in energy and momentum: w(po, |p|) ~ wo ©(po — E F ) 0(|p| — p F )- Here p F 
represents the nucleon Fermi momentum and E F = p F /2m N is the corresponding Fermi 
energy. By requiring this approximation to reproduce the observed flfov decay rates, we 
find wq ~ 4 MeV -1 , and by requiring that the present upper limit on the half life for 
neutrinoless decay [15] imply an upper limit for light neutrino masses [16] of \U 2 v m v \ < 2 
eV, 2 we find p F ~ 60 MeV (and so E F ~ 2 MeV). 

With these choices we obtain the following expression for the total PPo^ decay rate: 



r(/W- — ^ — \m , 



(15) 



with W given by an elementary integral: 



47T 3 



du 







p F — \l u 2 + m 2 tan 1 
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+ 



E- 



3m 7 - 



+ 



(E F < p F < m^). 



(mi <C E F <C p F ), 



(16) 



Motivated by the expression for the decay rate due to light neutrinos, it has become 
conventional to quote the experimental limit on f3f3$ v decay as an upper limit on the 
'effective' neutrino mass, defined by rrf^ = Yli^l^i- mentioned previously, the 
current experimental limit [15] implies an upper bound rrfj^ < 2 eV [16]. From the above 
formulae for T((3(3ou) we can obtain similar constraints on U e i and rrii for arbitrary neutrino 
masses. For the model under study here, we assume the v' e mass to be too small to 
contribute, and so there are three limiting cases to consider: 



We take here a bound which is twice as large as the usually-quoted limit, since we allow for an 
uncertainty in the nuclear matrix elements of a factor of 2. 
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(1) If the two sterile neutrinos are both light compared to a few MeV, then one would 
expect ml s = s\ (m n > + c% — M N >_ s^j as the approximate expression for m^, ff . However, 
in the model considered here this quantity vanishes, as may be seen from eq. (9). 
As a result we must work to sub-leading order in the sterile neutrino masses when 
approximating the integral in eq. (16), leading to 



m 



eff 



7TSi 



1 , 2E F 
- + log 



1 , 2E F 
- + log 

2 I Mat' 



< 2 eV. 



(17) 



(2) If both sterile neutrinos have large masses compared to p F ~ 60 MeV, then we instead 
find 



2 2 



M N ' + M N[ 



< 2 eV. 



(18) 



(3) Finally, if M N , + < p F and M N ,_ > p F: then 



M N , + 4 - 



2 2 



3M N , 



< 2 eV. 



(19) 



With these estimates we can determine the masses and couplings of the sterile neu- 
trinos which are consistent with the present non-observation of I3f3 0v . The constraints we 
obtain for U ei and m; in this way are plotted in Figs. (1). The area below the curves in 
this figure represents the allowed range. To see if a sterile neutrino will make observable 
contributions to l3f3o u decay, we must see if the range of values which lie close to those 
curves are consistent with other constraints. 



3.1) Radiative Corrections 



Our estimate of the (3(3q v decay rate in this section assumes a negligible contribution 
from the exchange of the very light neutrino mass eigenstate, v' e . This assumption requires 
some justification in parts of the parameter space which we consider here. We therefore 
briefly pause to provide this justification. 

The neglect of the light neutrino contribution to /3/3 u relies on our decision to choose 
this mass to be much lighter than 1 eV. At tree level this can always be accomplished by 
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choosing the parameters A e and n + to be sufficiently small. A naturalness problem can 
arise, however, if is chosen to be too small in comparison with the parameter 
which provides the lepton-number violating contribution to the heavy-neutrino masses. 
This is because a nonzero \i- generates, through loops, nonzero contributions to A e and 
, and if these loop-induced contributions are large enough, then the v' e exchange graph 
can only be neglected in (3(3o v if the loop-induced mass is cancelled by the tree-level term. 
A naturalness problem arises if the required cancellation becomes implausibly precise. 

More quantitatively, we imagine our model to be an effective theory which is obtained 
after some unknown physics above some scale A has been integrated out. We then use the 
renormalization group to run the couplings in this effective theory down from the scale A 
to the much lower energies that are relevant for (3(3 decay. In this way we can compute 
the contribution to the v' e mass which is produced by the mixing between the parameters 
A e and and fi- as they are run down from the scale A. We regard the theory to be 
natural if these contributions to the v' e mass are not larger than, say, 0(1 eV), and so do 
not need to be carefully cancelled by a 'bare' contribution to m v > at the scale A. 

The dominant graph to consider is that of Fig. (2), in which a SM Higgs scalar is 
emitted and absorbed by the light neutrino state. Its contribution to the light-neutrino 
mass is, in order of magnitude: 



For numerical purposes we take the logarithm in this expression to be unity. Requiring 
the rest of the result to be smaller than 0(1 eV) then produces the bound labelled NN 



We now turn to the exploration of the other constraints on this model. We consider 
in turn the limits coming from (i) supernova SN1987a, (ii) nucleosynthesis, and (Hi) 
laboratory limits. 




(20) 



in Fig. (3). 



4. Phenomenological Constraints 
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4.1) Supernova SN1987a 



The first constraint we consider is for low- mass sterile leptons ('low- mass' here means 
masses smaller than ~ 50 MeV) that mix significantly with v e . Any such particle can 
be produced in the core of a supernova, where temperatures are typically of order T SN ~ 
(60 — 70) MeV. This must be forbidden since otherwise the supernova would cool too much 
to be in agreement with the observations of SN1987a. The analysis we require is very 
similar to the case of right-handed neutrinos discussed in Ref. [17], where it was shown 
that these considerations lead to the following bounds 

3 x 10" 2 < |*7ei|, or else \U ei \ < 10" 5 . (21) 

The lower bound comes because for sufficiently strong mixing, the produced sterile neu- 
trinos get trapped in the supernova and so they do not provide a mechanism for cooling 
too quickly. The upper bound comes from the requirement that not too many sterile neu- 
trinos be radiated by the ordinary particles in the supernova. These bounds are displayed 
in Fig. (3) by the vertical and horizontal lines labelled SN. The vertical line is due to 
this bound being independent of the neutrino mass, provided only that this mass is much 
smaller than T SN . For the models of interest here, the upper bounds we obtain in this way 
for Mi and \U e i\ imply that the effective mass, m^, ff , which appears in (3(3q u decay can be 
at most rrf^ < 5 x 10~ 3 eV. This is well beyond the reach of present- and next-generation 
[18] /3/3-decay experiments. It is, however, worth noting that for M N > > 50 MeV, this 
bound is ineffective, and so SN1987a cannot rule out a significant contribution of sterile 
neutrinos to the (3(3qu decay rate, provided that these neutrinos are in this larger-mass 
regime. 

We therefore now turn to the constraints on N± which apply if the masses are 50 MeV 
or higher. 

4.2) Cosmological Constraints 

One of the major triumphs of the standard hot big-bang model of cosmology is its 
ability to naturally explain the primordial abundance of the light elements — 4 He, 7 Li, 
2 H and 3 He — using the standard model (SM) of electroweak interactions [19]. This, in 
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turn, implies stringent bounds on any new physics beyond the SM that involves weakly- 
interacting particles. In particular, in the model of present interest, the sterile neutrinos 
can upset the success of the nucleosynthesis discussion unless their masses and mixings are 
suitably constrained. 

The basic condition is to ensure that the energy density at the nucleosynthesis tem- 
perature, T BBN ~ 1 MeV, due to the sterile neutrinos is much less than that of an ordinary 
neutrino species. For sterile neutrinos that are much lighter than 1 MeV, this can be ar- 
ranged simply by having them decouple early enough for their energy density to be diluted 
by the reheating of ordinary matter, such as in the QCD phase transition. But keeping in 
mind the supernova constraint of the previous section, the sterile neutrinos in the model 
of interest here satisfy M 3> T BBN , and so even if they decouple sufficiently early, their 
relic energy density at nucleosynthesis will nevertheless dominate that of an ordinary neu- 
trino species unless their lifetime is shorter than 0.1 sec. Assuming decays through the 
charged-current weak interactions, this implies: 



Suppose, first, that the sterile neutrino is relativistic at the time that it decouples: 
T n > M N > . This is the case if 



Both of these conditions must be satisfied by both of the heavy particles, and are plotted 
as curves a and b in Fig. (3). BBN demands that any particle which lies to the left of 
curve b {i.e. decouples relativistically) must also lie to the right of curve a. Notice that 
these conditions together drive one, in the model of current interest, into a regime which is 
excluded by the laboratory bounds considered in the following section, provided one asks 
for a I3{3q v rate close to observability at the same time, and so the relativistic-decoupling 
scenario is not relevant in our case. 

Next suppose that the sterile neutrinos decouple nonrelativistically, T D < M N > ± , and 
so lie to the right of curve b in Fig. (3). In this case their number density is exponentially 




(22) 



± 




(23) 
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suppressed by the Boltzmann factor. This dilution then makes them decouple earlier 
than one would otherwise expect without the Boltzmann suppression. The decoupling 
temperature, T D , is therefore given by 



M 2 N± G 2 F \U ei 



(M JV± T D ) 3 / 2 e~ Mjv±//TD < gV 




(24) 



TV 



M pl 



In the case of present interest the sterile neutrinos will always decouple at or above a 
temperature of ~ lOOMeV, i.e. lie to the right or below curve c in Fig. (3), provided flflov 
occurs at or close to an observable level. This in turn means that their energy density at 
nucleosynthesis will only be sufficiently small if they, again, decay fast enough, eq. (22). 

To summarize, having observable /3/3 0l , in this model implies that the heavy sterile 
neutrinos decouple only after having become non relativistic, but still early enough that 
conflict with standard Big-Bang Nucleosynthesis can be avoided, provided they decay fast 
enough. 

4.3) Laboratory Limits 

A wide range of experiments constrain the properties of isosinglet heavy leptons [2] , [3] , 
[4], [5], [6], [7]. They can be classified into two main categories, according to whether the 
bound is obtained from precision measurements on the Z resonance, or from experiments 
at lower energies. 

There are two main types of low-energy experiments which limit the properties of 
sterile neutrinos that mix with the electron neutrino, such as for the model considered 
here. One type obtains its bound from the decay rate and the electron spectrum of the 
two-body decay of kaons and pions at rest. For example, for neutrinos in the mass range 
between 1 and 100 MeV, the measured r(/T — > ev)/T(n — > [iv) rate provides a mass- 
dependent bound on \U e {\. For a 1 MeV neutrino the bound is \U e i\ 2 < 10~ 3 at the 90% 
C.L., whereas the respective bounds for a 10 MeV and a 50 MeV neutrino are 10 -5 and 
5 x 10 -7 [20]. Similar searches for a nonstandard component to K — > ev [21] extend this 
limit up to sterile-neutrino masses of 350 MeV. 

Even stronger limits can be obtained from beam-dump experiments provided that the 
heavy neutral leptons can decay appreciably through their charged- current interactions. 
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These can constrain sterile- neutrino masses up to ~ 2 GeV, with a sensitivity of \U ei \ 2 < 
10~ 7 . For masses below ~ 0.5 GeV the bound becomes as good as \U ei \ 2 ~ 0.5 x 10~ 9 
[22]. (For a more detailed discussion see [3]). 

For neutrinos with masses that are more than a few GeV, measurements at the Z 
pole provide the strongest limits. For M N > < M z , the best bounds come from the nonob- 
servation of the decay of a Z into a sterile and a standard neutrino, Z — > Nu — > W*ev, 
with the subsequent decay of the sterile neutrino through a virtual boson, W*. The bound 
obtained in this way is \U ei \ 2 < 7 x 10 -5 [2], [4], [7]. 

The above bounds do not apply if the heavy singlet neutrino is heavier than M z . In 
this case there are two types of bounds to consider, which arise due to the reduction of 
the couplings of the ordinary neutrinos to the W and Z bosons due to their admixture 
with the new sterile neutrinos. The reduction in the effective couplings to the Z, result 
in a reduction of the Z : s invisible width. This leads to the bound: \U ei \ 2 < 2.7 x 10~ 2 
for isosinglet masses above 90 GeV [6], [7]. A stronger bound arises however from the 
reduction of the W couplings, which potentially show up as a failure of lepton universality 
in low-energy weak decays, as well as affecting precision electroweak measurements through 
their influence on the experimental value of Fermi's constant, G F , that is inferred from 
muon decay. This implies the bound : \U e i\ 2 < 5.6 x 10~ 3 (2a) for isosinglet masses above 
90 GeV [6]. 

All the phenomenological constraints discussed in this section have been summarized 
in Fig. (3). When contrasted with the masses and mixing angles required for a /3/3 0u 
signal close to observability in the model discussed here (Figs. (1)) they yield the region 
of parameter space depicted in Fig. (4) . 

5. Conclusions 

Our purpose here has been to determine which kinds of heavy sterile neutrinos can 
contribute appreciably to PPq u decays, and to explore the constraints which such particles 
must satisfy due to present laboratory and astrophysical information. Part of our moti- 
vation for doing so has been to provide an example of a theory in which this decay can 
proceed without requiring the existence of a light neutrino with a mass in the vicinity of 
1 eV. As an existence proof for theories of this type, we display here a model which does 
so. It does so by producing observable Pfiov purely through the exchange of a sterile neu- 
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trino having a mass in the GeV range. This runs contrary to the usual expectation that 
the observation of f3f3 0l/ must indicate the existence of a majorana mass for the electron 
neutrino in the eV range. 

The model we consider is reasonably simple, supplementing the standard model only 
by two new left-handed neutrino states. We find that requiring the model to be consistent 
with all astrophysical and laboratory limits, as well as with an observable (3(3q u signal, 
constrains the couplings and masses of the new neutrinos to lie in a limited region of 
parameter space. The mass range that is favoured by these bounds, as well as naturalness 
considerations, is 1 - 10 GeV. For these masses, the couplings that are required to produce 
an observable PPov signal are roughly \U e i\ 2 ~ 10~ 5 . Such parameters place such a sterile 
neutrino close to the current limits of detection at LEP, where they can be searched for 
through the decay Z — > iW e , with the subsequent charged-current decay of the sterile 
neutrino, N, into quarks and leptons. This shows how LEP results can be used to help 
diagnose the implications of a potential /3/3 0l , signal, and illustrates the rich interplay that 
is possible between low- and high-energy experiments. 
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7. Figure Captions 



This figure plots the region in the U ei — mi plane (Fig. la, (lb) for N' + , (N'_ respec- 
tively)) that is allowed by current (3(3q u experiments for various different values of 
the parameter e as defined in eq. (8). The area below the corresponding lines is the 
allowed region. The limiting case of very large e appears as the solid line labelled 
"e = oo" , which also corresponds to the contribution of a single sterile neutrino that 
mixes with v e . 

The Feynman graph which dominates the renormalization-group mixing of with 
A e and fx + . 

This figure summarizes the phenomeno logical constraints as discussed in section 4. 
The region above the solid line is ruled out by the various laboratory bounds. The 
region below and to the left of the dashed line labelled SN is excluded by the obser- 
vations of SN1987a. The dashed line labelled NN depicts the naturalness bound as 
discussed in section 3 (Eq. (20)) for e = 0.1. In the region to the right and above 
this line fine tuning is required. Finally the dash-dotted curves labelled a, b and c 
represent the nucleosynthesis bounds. The lifetime of a sterile neutrino is less than 0.1 
sec to the right of curve a. A particle decouples after having become non relativistic 
in the region to the right of line 6, but it will nevertheless decouple at or above a 
temperature of ~ 100 MeV to the right and below line c. Line a thus represents the 
only relevant nucleosynthesis bound in our case, and the region to the right of line a 
and below the solid line is allowed. 

We plot here the region of the parameter space in the e — M N i plane which is obtained 
by requiring masses and mixing angles which yield a flflov signal close to observability, 
as depicted by the various lines in Fig. (1), together with the various phenomenological 
bounds displayed in Fig. (3). The allowed area is marked by shading, and extends 
upwards beyond the region depicted in the figure towards higher values of e without 
changing the mass range. The darker area represents the part of the parameter space 
in which the smallness of the v' e mass is explained in a way which is technically 
natural in the sense explained in section (3.1) of the text. In the lightly-shaded region 
a finetuning of A e and fi+ is required in order to keep the v' e mass below 1 eV. Notice 
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that, since M N >_ > M N > , the analogous figure in the e — M N >_ plane would be shifted 
to the right according to eq. (9). 
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This figure "figl-l.png" is available in "png" format from: 



http://arXiv.org/ps/hep-ph/9408367v2 
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